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Abstract: The Axelrod model is a spatial stochastic model for the dynamics of cultures 
which, similarly to the voter model, includes social influence, but differs from the latter by also 
accounting for another social factor called homophily, the tendency to interact more frequently 
with individuals who are more similar. Each individual is characterized by its opinions about a 
finite number of cultural features, each of which can assume the same finite number of states. 
Pairs of adjacent individuals interact at a rate equal to the fraction of features they have 
in common, thus modeling homophily, which results in the interacting pair having one more 
cultural feature in common, thus modeling social influence. It has been conjectured based on 
numerical simulations that the one-dimensional Axelrod model clusters when the number of 
features exceeds the number of states per feature. In this article, we prove this conjecture for 
the two-state model with an arbitrary number of features. 



Oh '. 1- Introduction 

Oh 

r £^ . There has been in the past decade a rapidly growing interest in agent-based models in an attempt 

to understand the long-term behavior of complex social systems. These models are characterized 
by heuristic rules that govern the outcome of an interaction between two agents, and a graphical 
structure, modeling either physical space or a social network, that encodes the pairs of agents that 
may interact due to, e.g., geographical proximity or friendship. The main objective of research in 
this field is to deduce the macroscopic behavior that emerges from the microscopic rules, which 
also depends on the structure of the network of interactions. The mathematical term for agent- 
based models is interacting particle systems though, as scientific fields, the former involves numerical 
simulations whereas the latter is based on rigorous mathematical analyses. While there is a common 

t^- ■ effort from sociologists, economists, psychologists, and statistical physicists to understand such 

models, interacting particle systems of interest in social sciences have been so far essentially ignored 
by mathematicians, with the notable exception of the voter model. This paper is motivated by this 
lack of analytical results and continues the study initiated in [5] for one of the most popular models 
of social dynamics: the Axelrod model [2]. The effort to collect analytical results is mainly justified 

rS \ by the fact that stochastic spatial simulations are generally difficult to interpret. This is especially 

true for the Axelrod model which, in contrast with the voter model, has a number of absorbing 
states that grows exponentially with the size of the network, and for which simulations of the finite 
system can freeze in atypical configurations, thus exhibiting behaviors which are not symptomatic 
of the long-term behavior of their infinite counterpart. 

The Axelrod model [2] has been proposed by political scientist Robert Axelrod as a stochastic 
model for the dissemination of culture. The heuristic microscopic rules include two important social 
factors: homophily, which is the tendency of individuals to interact more frequently with individuals 
who are more similar, and social influence, which is the tendency of individuals to become more 
similar when they interact. Note that the voter model [3, 4] accounts for the latter but not for 
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the former: individuals are characterized by one of two competing opinions which they update at 
a constant rate by mimicking one of their neighbors chosen uniformly at random. In particular, 
any two individuals in the voter model either totally agree or totally disagree, which prevents 
homophily from being incorporated in the model. In contrast, individuals in the Axelrod model 
are characterized by a vector, also called a culture, that consists of F coordinates, called cultural 
features, each of which assumes one of q possible states. Homophily can thus be naturally modeled 
in terms of a certain cultural distance between two individuals: pairs of neighbors interact at a 
rate equal to the fraction of features they have in common. Social influence is then modeled as 
follows: each time two individuals interact, one of the cultural features for which the interacting 
pair disagrees (if any) is chosen uniformly at random, and the state of one of both individuals is 
set equal to the state of the other individual for this cultural feature. More formally, the Axelrod 
model on the infinite one-dimensional lattice, which is the network of interactions considered in this 
paper, is the continuous-time Markov chain whose state space consists of all spatial configurations 

r]:Z — > {l,2,...,q} F . 

To describe the dynamics of the Axelrod model, we let 



1 F 
F(x,y) = -^ 1{t](x, i) =r](y,i)}, 



where rj(x, i) refers to the ith coordinate of the vector r)(x), denote the fraction of cultural features 
vertex x and vertex y have in common. In addition, we introduce the operator <r x ,y,i on the set of 
spatial configurations defined by 

i \t -\ iv(y,i) tf z = x and j = i . 

{o- x ,y,iV){z,3) = < , ., ,, for x,y G Z and i G {1,2,... ,q\. 

I rj{z,j) otherwise 

In other words, configuration cr x ,y,i V is obtained from configuration rj by setting the ith feature of 
the individual at vertex x equal to the ith feature of the individual at vertex y and leaving the 
state of all the other features in the system unchanged. The dynamics of the Axelrod model are 
then described by the Markov generator L defined on the set of cylinder functions by 



L fw = E E ^ 



F(x,y) 



2F 

\x-y\ = l «=1 



1-F(x,y) 



l{y(x,i) i- V(y,i)} [f(vx,y,iV) - fin)]- 



Note that the expression of the Markov generator indicates that the conditional rate at which the 
ith feature of vertex x is set equal to the ith feature of vertex y given that these two vertices are 
nearest neighbors that disagree on their ith feature can be written as 



_L_ f F(x,y) 
2F [l-F(x,y) 



1 1 

F(x,y) x — ; — x - 

1 ,yj F(l-F(x,y)) 2 



which, as required, equals the fraction of features both vertices have in common, which is the rate 
at which the vertices interact, times the reciprocal of the number of features for which both vertices 
disagree, which is the probability that any of these features is the one chosen to be updated, times 
the probability one half that vertex x rather than vertex y is chosen to be updated. 
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The main question about the Axelrod model is whether or not the population converges to a 
consensus when starting from a random configuration. For simplicity, we assume that the initial 
cultures at different sites are independent and identically distributed, and that at a given site, 
each of the q possible initial cultures appears with the same probability. The term "consensus" is 
defined mathematically in terms of clustering of the infinite system: the model is said to cluster if 

lim P (i]t(x, i) = T] t (y, i)) = 1 for all j;,t/6Z and i G {1, 2, . . . , g} 

t— >oo 

and is said to coexist otherwise. The dichotomy between clustering and coexistence for the finite 
model is unclear since, as mentioned above, the finite system can hit an absorbing state in which 
different cultures are present even though its infinite counterpart clusters. In order to characterize 
the transition between the two regimes, Vilone et al [6] considered the random variable s max which 
refers to the length of the largest interval in which all individuals share the same culture in the 
absorbing state hit by the finite system, and distinguished between the two regimes depending on 
whether the expected value of this random variable scales like the population size or is uniformly 
bounded. Denoting the population size by N, their spatial simulations suggest that 

limTv^oo E (iV _1 s max ) > for F > q and limTv^oo E (s max ) < oo for F < q 

so we conjecture clustering when F > q and coexistence when F < q for the one-dimensional 
infinite system. The mathematical analysis of the Axelrod model initiated in [5] strongly suggests 
the coexistence part of the conjecture for a certain subset of the parameter region. More precisely, 
letting n c denote the number of cultural domains in the absorbing state hit by the finite system, it 
is proved based on duality-like techniques and a coupling with a simple urn problem that 

liniAr^oo E (A r ~ 1 n c ) > for F < c x q where c « 0.567 satisfies c = e~ c . 

It is also proved that the infinite system clusters in the critical case F = q = 2. In this paper, we 
extend this result to all values of the number of features when the number of states per feature is 
again equal to two, which proves the clustering part of the conjecture stated above when q = 2. 
More precisely, we have the following theorem. 

Theorem 1 - The F-feature 2-state Axelrod model on Z clusters, starting from the random initial 
configuration in which the cultures at different sites are independent and identically distributed 
and at a given site, each of the 2 F possible initial cultures appears with the same probability. 

The proof when F = q = 2 is carried out in [5] based on duality techniques for the voter model 
through the existence of a natural coupling between the two-feature two-state Axelrod model and 
the voter model obtained by identifying cultures with no feature in common. This coupling, however, 
fails for any other values of the parameters, so a different approach is needed to extend the result to 
a larger number of features. The first step is to construct a coupling between the two-state Axelrod 
model and a certain collection of non-independent systems of annihilating symmetric random walks 
that keep track of the disagreements between nearest neighbors. Clustering of the Axelrod model 
is equivalent to extinction of these systems of annihilating random walks. The proof of the latter 
is inspired by a symmetry argument introduced by Adelman [1] which is combined with certain 
parity properties of the collection of non-independent systems of random walks. 
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Figure 1 . Pictures of the system of coalescing random walks obtained from a realization of the two-state Axelrod 
model on the one-dimensional torus with 600 vertices with F — 3 features. The left picture gives a superposition of 
all three levels of random walks while the first level only is extracted in the right picture. 



2. Systems of annihilating random walks 

In this section, we represent the Axelrod model on Z by a particle system that keeps track of 
the interfaces between cultural domains, thus looking at the disagreements along the edges of the 
graph rather than the actual cultures on the vertices. This approach is motivated by the fact that 
consensus in the Axelrod model is equivalent to the extinction of its interfaces. In order to obtain 
a well-defined Markov process, it is necessary to keep track of the features for which neighbors 
disagree rather than simply the number of these features. Therefore, we think of each edge of the 
graph as having F levels, and place a particle on edge e = (x, y) at level i if and only if vertex x 
and vertex y disagree on their ith feature. That is, we define 



&(u,t) = l{Tfr(u-l/2, i)^n t (u + 1/2, «)} for all u£ 



+ 1/2 



and place a particle at site u € B at level i whenever £t(u, i) = 1. To study the dynamics of this 
system and the rates at which particles jump, it will be useful to also keep track of the number of 
particles per site so we introduce 



Ct(u) = ^ ^( n '*) for all n e 



and call u a j-site whenever it has a total of j particles: Ct(u) — j- To understand the evolution, 
we first observe that when a particle jumps, it moves right or left with equal probability, unless 
another particle already occupies the site on which the particle tries to jump in which case the 
particles annihilate each other. Thus, these processes induce a collection of F non-independent 
systems of annihilating symmetric random walks. The symmetry is due to the fact that, when two 
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neighbors interact, each of them is equally likely to be the one chosen to be updated. Also, the 
reason why a collision between two particles results in an annihilation is that when two individuals 
disagree with a third one on a given feature, these two individuals must agree on this feature, which 
happens when q = 2. (Note however that for larger values of q collisions between particles would 
result in either an annihilating event or a coalescing event depending on the configuration of the 
underlying Axelrod model.) Even though the evolution of the particle system at a single level is 
somewhat reminiscent of the evolution of the interfaces of the one-dimensional voter model, it is 
in fact much more complicated due to the presence of strong dependencies among the different 
levels. These dependencies result from the inclusion of homophily in the model, which implies that 
particles jump at varying rates. More precisely, since two adjacent vertices that disagree on exactly 
j of their features interact at rate 1 — j/F, the fraction of features they share, and the site between 
these vertices is a j'-site, given that u is a j'-site each particle at site u jumps at rate 

r(j) = --— for j ^ 0. 

This indicates that the motion of the particles is slowed down at sites that contain more particles, 
with the dynamics being frozen at sites with F particles. In the following sections, we will call 
particles frozen or active depending on whether these particles are located at an F-site or not, 
respectively. We refer the reader to Figure 1 for simulation pictures of the systems of annihilating 
random walks when F = 3. 



3. Showing that the process can not become frozen 

In this section, we prove that no site can remain an F-site forever, which is the key to proving 
consensus in the two-state Axelrod model. From the point of view of the particle system described 
in the previous section, this means that if some site u is completely filled with particles at some 
time, then eventually another particle will jump onto the site u, annihilating one of the frozen 
particles on the site and making the other F — 1 particles active. 

Proposition 2 - Assume that site u € ID is a F-site at time t. Then, 

T := inf{s > t : ( s (u) ^ F} < oo a.s. 

It will be useful in the proof of the proposition to construct graphically the particle system described 
in the previous section using the following collections of independent Poisson processes and random 
variables: for each pair of site and feature (v, i) € B x {1, 2, . . . , F}, 

• we let (N v j(t) : t > 0) be a rate one Poisson process, 

• we denote by T Vt i(n) its nth arrival time: T V i(n) = inf {t : N V j(t) = n}, 

• we let (B v j(n) : n > 1) be a collection of independent Bernoulli variables with 

P(J3 W)i (n) = +l) = P (B Vti (n) = -1) = 1/2, 

• and we let (U Vt i(n) : n > 1) be a collection of independent Uniform (0, 1) random variables. 
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The system of annihilating random walks is constructed as follows. At time t = T Vy i(n), we draw 
an arrow labeled i from site v to site v + B v ^{n) to indicate that if 



£t-(v,i) = 1 and U Vji {n) < r(( t -(v)) 



1 1 



Ct-(v) F 



then the particle at site v at level i jumps to site v + B V j(n). 

The above construction can be extended naturally to any subgraph G of the lattice by using the 
same collections of independent Poisson processes and random variables and killing all the particles 
that jump onto a site v which is not the center of an edge of the graph. Consider now the case in 
which the graph G is the one induced by the vertex set N. Suppose the initial configuration is such 
that the left-most edge has F particles, one at every level, while every level of every other edge 
independently has a particle at time zero with probability 1/2. Let 

p := P (the left-most edge has F particles at all times). 

That is, p is the probability that no particle ever tries to jump onto one of the particles on the 
left-most edge. We will later see that p = 0. 

Returning now to the setting of Proposition 2, fix a site «eD and a time t > 0, and suppose 
that particles at that space time point are frozen: Ct(u) = F. We will consider only the sites to the 
right of u and show that eventually some particle from the right must jump onto u, unless the site 
u has already been hit from the left. Lemma 5 below will allow us to break the process into stages, 
and give a lower bound for the probability that a particle at u is annihilated at each stage. We 
start by proving Lemma 4 which is a key preliminary result. 

Definition 3 - The interval {u, u + 1, . . . , v} CD is said to be active at time t if the numbers of 
particles it contains at two different levels differ in their parity, i.e., 

V V 

^2 &( w >*) ^ 5Z &( W >J) mod 2 for some i^ j. 

w=u w=u 

Lemma 4 - Assume that {n, . . . , v} is active at time t and Ct(u) = Ct(v) = F. Then, 

T := inf{s > t : ( s (u) + ( s {v) / 2F} < oo a.s. 

Proof. Seeking a contradiction, we assume that Cs{ u ) = Cs( v ) = F for all s > t. Under this 
assumption, the parity of the number of particles at each level between site u and site v is preserved 
since particles annihilate by pair. In particular, {u, . . . , v} is active at every later time, which implies 
that this interval contains at least one site which is neither a 0-site nor an F-site, and thus at least 
one active particle, at every later time. Since this particle jumps at a positive rate, it must hit one 
of the boundaries u or v in an almost surely finite time, which leads to a contradiction. □ 

Lemma 5 - There exists a sequence of times t = to < t\ < ti < ■ ■ ■ < oo such that, if A^ denotes 
the event that at some time s € (tk— i,tk], a particle at site u is annihilated, then 

P[A k \ f]A%) > | for all k > 0. (1) 
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Proof. The proof relies in part on delicate symmetry arguments and is modeled after the con- 
struction of Adelman [1]. We refer the reader to Figure 2 for a picture of this construction in our 
context. We must analyze in detail the process at each stage. For w € D, w > u, let 

F w {t) = a ((N v>i (s) :0<s<t), (B„,<(n) : 1 < n < N v>i {t)), 

(U V)i (n) : 1 < n < N v>i (t)) : u < v < w, I < i < q) 

be the <r-field generated by the graphical representation of the systems of random walks over the 
spatial interval {u, u + 1, . . . ,w} through time t. We will prove (1) by induction on k. As part of 
the induction hypothesis, we will assume that at the beginning of the fcth stage of the construction, 
we have a site v k such that the following two conditions hold: 

HI. N Vh>i (t fc _i) = N Vk+1>i (t fc _i) = N Vk+2 ,i (tfc-i) = for i = 1, 2, . . . , F, 
H2. A 1 ,A 2 ,...,A k - 1 er Vk (t k -i). 

For k = 1, condition H2 is trivial while HI can be satisfied by choosing to = t and 

V! = mm{v > u : N V)i (t) = N v+1>i (t) = N v+2 ,i(t) = for i = 1, 2, . . . , F}. 
To prove (1) by induction, we will show that 

P(A k \F Vk (t k ^)) > | (2) 

and that t k and v k +i can be chosen to satisfy conditions HI and H2. 

Condition HI above implies that, in the graphical representation, there is no arrow starting at 
either site v k or site v k + 1 by time t k -\ from which it follows that 

• particles starting to the right of v k do not reach v k by time tk-i, and 

• particles starting in {u, u + 1, . . . , v k } do not reach v k + 1 until after time t k -i- 

Note that the assumption N Vk+2 ,i (tk—i) = is not necessary at this stage of the proof but it will 
be useful later in the proof to obtain a lower bound of the probability that a certain interval is 
active. We will partition the region to the right of site v k + 2 into intervals of equal length, each 
containing a total of v k — u + 3 sites. Specifically, for each positive integer j, we let 

Bj,k = {J v k ~ (j ~ 1)« + 3j, . . . , (j + l)v k - ju + 3j + 2}. 

Also, we define the interval B ok = {u, u + 1, . . . , v k + 2} as depicted in Figure 2. Now let J k be the 
smallest positive integer j such that the following two conditions hold: 

• we have N Wy i(t k -.i) = for all w € Bj k and i = 1, 2, . . . , F, and 

• we have either 

translation : ^(jv k - (j - l)u + 3j + m,i) = & fc _i (« + m , i) 
or reflection : £ Q ((j + l)v k - ju + 3j + 2 - m, i) = Ct fc _i(w + m,i) 

for m = 0, 1, . . . , v k — u + 1 and i = 1,2, . . . ,F. 
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The second requirement ensures that either the configuration in Bj k ^ excluding its right-most site at 
time zero is the same as the configuration between u and v k + l at time tk-i, or else the configuration 
in Bj k ^ excluding its left-most site at time zero is the mirror image of the configuration between 
u and vj- + 1 at time tfc-i- In the latter case, we say that reflection occurs, an event which, due to 
obvious symmetry, has probability 1/2. 

We now consider the process conditioned on the event J k = j as well as the event that reflection 
occurs. We also, for now, truncate the process at the right edge of Bj t k, so the right-most site is 

rk ■= (j + 1) Vk - ju + 3j + 2 

and just consider the evolution of the process between the sites u and r k - Let 

4 := j v k - {j - 1) u + 3j 

be the left-most site in Bj k . We make the following observations: 

• Because no particles in the interval Bj k can jump before time tk—i, the left-most and the 
right-most vt — u + 1 sites are mirror images of each other at time tk-i, i.e., 

ft*.! (u + m, i) = &*_i(nfc - Tn,i) for m = 0, 1, . . . ,v k - u+ 1 and alii 

• Because particles at the sites Vk, v k + 1, v k + 2, £)., and £k + 1 do not jump before time tk-i, 
conditional on the configuration at the left-most and the right-most v^ — u + 1 sites at time 
tfc_i, the sites in between evolve independently from time until time tfc-i- The law of their 
evolution is the same as the law of the original process, modified so that particles at the first 
two sites and the last two sites are not permitted to jump and conditioned on the blocks Bi^ 
for i = 1, . . . , j — 1 failing to satisfy the conditions of either translation or reflection. As a 
result of the way </& is defined, the law of this process is therefore the same as the law of its 
mirror image. 

• Because particles at site Vk + 1 and site 4 + 1 have no opportunity to interact with other 
particles before time tk—i, on these sites at each level there is a particle independently with 
probability 1/2. Consequently, the probability that the number of particles at two given levels 
of these two sites have the same parity is equal to 1/2. 

By the last observation above, if a reflection occurs, then the interval from u to r/, is active in 
the sense of Definition 3 with probability at least 1/2. Furthermore, by the symmetry noted in 
the second observation above, the first change either to site u or site r^ happens to site u with 
probability 1/2. Hence, if Dk denotes the event that eventually there is a change to site u, then 



P{D k \F Vk {t k ^)) > 



IV 



To complete the construction, we let £& = tk-i in case a reflection does not occur or the interval 
from u to rfc is not active. Otherwise, we define time t k to be the first time at which there is a 
change either to site u or to site r/% in this system. 

We now return to the case in which only the sites to the left of u are discarded, which means 
there is a possibility that some particle could jump onto r k from the right. In this case, conditional 
on the event of a reflection, which implies that the particles at r k are frozen until at least t k , and 



u 



O 

cx: 

CO 
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v k + 2 l k 



r k 



B 0,k 


B l,k 


J k = 2 
mirror of Bq k 










J k-l = 3 


-* •— • • ►- 



f k 



tk-1 



t u*-i + 2 t 

frozen no jump before t k -i 

Figure 2. Picture related to the proof of Lemma 5 



no jump before t k - 



frozen 



thinking of site r k as the left-most edge in the definition of p, we see that there is a probability 
larger than p that no particle will jump from r k + 1 to r& before time t k - We deduce that there is 
probability at least p/8 that site u will change by time t k , which gives (2). Finally, we let 



Vk+l 



min {j > r k : N V)i (t k ) = N v+lji (t k ) = N v+2ji (t k ) = for i = 1, 2, . . . , F}. 



Because no particle to the right of v k +i can reach v k by time t k , we have A k € F Vk+1 {tk), which 
completes the proof of the lemma. □ 

Proof of Proposition 2. Consider the system in which sites to the left of u are removed. We 
proceed by contradiction by assuming that there is a positive probability that the site u never 
changes after time t. We will show that this implies that p > 0. Thus, (1) will imply that 



p(QV) 



1 



which means that with probability one, eventually some particle at site u will be annihilated. This 
contradiction will imply the result. 

In the case t = 0, the probability that the site u never changes after time t is exactly p, by 
definition. We still need to show that if site u becomes an F-site at t > 0, then having a positive 
probability that the site never changes after time t would still imply that p > 0. The strategy will 
be to argue that with positive probability, up to time t the process behaves within a finite region 
in such a way that makes it possible to reduce to the t = case. 

Because, with probability one, there are infinitely many sites v > u such that 

N v ,i(t) = N v+1>i (t) = fort = 1,2,...,*; 

there exists a site v > u such that 



p (C s ( u ) = F for all s > t and N vi (t) = N v+li (t) = for i = 1, 2, . . . , F) > 0. 
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The event N v ^(t) = N v +n(t) = for i = 1,2, . . . ,F implies that there is no arrow connecting site 
v and site v + 1 by time t. Therefore conditional on this event, the evolution of the process on the 
interval {u, u + 1, . . . ,v} is independent of its evolution on {v + 1, v + 2, . . . } up to time t. 

Because there are only finitely many possible configurations for the sites u + 1, . . . ,v — 1, there 
must exist numbers c w j £ {0, 1} for u < w < v and i = 1,2, . . . F such that 

P (Cs(u) = F for all s > t, £,t(w, i) = c w ^ for u < w < v and i = 1,2, . . . , F, 

and N v ,i(t) = N v+lii (t) = for i = 1, 2, . . . ,F) > 0. 

Let Cu y i = 1 for i = 1,2, ... ,F. Clearly there is a positive probability that 

£,o(w, i) = c Wt i for w = u, u + 1, . . . , v and i = 1,2, . . . , F and 
N w> i(t) = for w = u, u + 1, . . . , v and i = 1,2, . . . ,F 

so it follows from the conditional independence noted above that 

P {( s (u) = F for all s > t, £,o(w, i) = c w ^ for u < w < v and i = 1,2, . . . , F, 

and N Wt i(t) = for w = u, u + 1, . . . , v + 1) > 0. 

However, the probability on the left-hand side is at most 

p (Q( u ) = F for all s > and £ Q (u, i) = 1 for i = 1, 2, . . . , F), 

which in turn is at most p. It follows that p > 0. □ 

4. Extinction of the particles 

In this section, we prove almost sure extinction of the systems of annihilating random walks. That 
is, we show that 

lim P(( t (u)=0) = 1 for all USB. (3) 

t— >oo 

We use different strategies to deal with active particles and frozen particles. We start by proving 
extinction of the active particles since it is one of the keys to showing extinction of the frozen 
particles, but the main ingredient to prove the latter is the result of Proposition 2. To see that (3) 
indeed implies Theorem 1, we observe that for all x, y £ Z with x < y, we have 

lim P (rj t (x) ^ r] t (y)) < lim P (r) t (x + j - 1) ^ rj t (x + j) for some j = 1,2, ... ,y - x) 

t— >oo t— >oo 

< lim P ((t(x + j — 1/2) / for some j = 1, 2, . . . , y — x) 

t— >oo 

y-x 



< lim V P (Ct(x + j - 1/2) ^ 0) = 0. 



t— )-oo 

i=i 



It remains to prove almost sure extinction of the systems of annihilating random walks, which is 
done in Lemma 7 for the active particles and in Lemma 8 for the frozen particles. Note that Lemmas 
7 and 8 immediately imply (3). 
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Lemma 6 - The limits 

lim E(Ct(u) l{(t(u) + F}) and lim E {Q t {u) l{C, t (u) = F}) 

t— >oo t— >oo 

exist and do not depend on the choice of u € ID. 

Proof. First of all, since the initial configuration as well as the evolution rules of the process are 
translation invariant in space, the probability that u has i particles at time t does not depend on 
the choice of a specific site «6B, i.e. 

p (Q( u ) =i) = P {C, t {v) = i) for all u, v € B. 

This implies that the limits superior and limits inferior of the expected values in the statements do 
not depend on u. In addition, since particles can only annihilate, the expected number of particles 
per site is a nonincreasing function of time so it has a limit as t — > oo and we write 

limsup E(( t {u)) = liminf E (Ct{u)) = lim E (£ t {u)) = L. (4) 

Now, seeking a contradiction, we assume that 

L_ = liminf E(( t (u) l{( t (u) ^ F}) < limsup E (( t (u) l{( t (u) ^ F}) = L+ (5) 

and let e > such that 5e = L + — L_. In view of (4) there exists to < oo such that 

L-e < E (Ct(«)) < L + e for all t > t . (6) 

Also, in view of (5) there exists an infinite sequence of times 

to < Si < t\ < S2 < t2 <■■■< Sj < tj <■■■< CO 

such that for all integers j > 0, we have 

E(( Sj (u) l{( Sj (u) ^F}) < L_ +e and E(( tj (u) l{Ct» + F}) > L+ - e. (7) 

It directly follows from (6) and (7) that 

E(( t] (u) l{Ct» = F}) - E(C..(u) l{C Sj (u) = F}) 

= E{Q tj {u)) - E(C Sj (u)) + E{C Sj {u)l{Q tj {u)^F}) - E (( tj («) 1{C S , («) ^ F}) 

< (L + e) - (L-e) + (L_ + e) - (L+ - e) = L_ - - L+ + 4e = -e 

for all j > from which we deduce that 

P(Ct j (u)=F) - P(C.(u)=F) 

(8 
= F- 1 E(Q tj (u)l{C t] (u) = F}) - F- 1 E(C Sj (u)l{( Sj (u) = F}) < -eF~\ 

Now, since each time an F-site becomes an (F — l)-site, there are two particles that annihilate each 
other, the inequality in (8) also implies that 

E{Q t .{u)) - E(CsAu)) < 2xP(CtM = F) -- 2xP(C Sl (u) = F) < -<leF~ x . 
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In particular, using once more that the expected number of particles per site is a nonincreasing 
function of time, and applying the previous inequality F times, we obtain 

F 

E(( tF (u)) ■- E(( Sl (u)) < ]T (£(&») -" E(( Sj (u))) < -2e, 

i=i 

which contradicts (6) above. In particular, (5) does not hold so the expected number of active 
particles per site has a limit as t — > oo. To prove that the second limit in the statement also exists, 
we simply observe that 

E(Ct(u) l{C t (u)=F}) = E(( t (u)) ■- E(( t (u) 1{0 < ( t (u) < F}) 

and invoke the fact that both expected values on the right-hand side have a limit as times goes to 
infinity. This completes the proof. □ 

Lemma 7 - There is extinction of the active particles, i.e., 

lim P (0 < ( t (u) < F) = for all u € D. 

t— >oo 

Proof. Note first that, by Lemma 6, the expected number of active particles per site does not 
depend on the choice of a specific site and has a limit as time goes to infinity. We assume by 
contradiction that this limit is positive, which implies that 

E (Ct(n) 1{0 < ( t (u) < F}) > M > for all t > 0. (9) 

We say that there is an event at site u at time t if one of the following occurs. 

1. Annihilating event: &(«) + (t(u + 1) = (t-(u) + ( t -(u + 1) - 2. 

2. Freezing event: Ct-(u) = F — 1 and C^^) = E. 

Denote by At(u) and Ft(u) the number of events at site u by time t, and observe that the joint 
distribution of these two random variables does not depend on the specific choice of the site u. This 
follows by again invoking the translation invariance in space of the initial configuration and the 
evolution rules of the process. Since active particles evolve according to symmetric random walks 
run at a positive rate, they are doomed to either annihilation or freezing. Hence, (9) implies that 

lim E(A t {u)) + lim E {F t {u)) = oo. (10) 

t— >oo t— >oo 

Now, we observe that, since the initial number of particles per site is bounded by F and each 
annihilating event removes two particles from the system, 

E(A t {u)) < F/2 < oo foralU>0. (11) 

In addition, two consecutive freezing events at site u must be separated by an annihilating event 
at either site u or site u — 1 which gives the upper bound 

E(F t (u)) < l + E{A t {u)) < l + F/2 < oo for all t > 0. (12) 

Since the combination of (11) and (12) contradicts (10), the lemma follows. □ 
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Lemma 8 - There is extinction of the frozen particles, i.e., 

lim P (Ct{u) = F) = for all ueB. 

t—>oo 

Proof. Note as previously that, by Lemma 6, the limit to be estimated exists and does not depend 
on the specific choice of site u£D. Let e > 0. By Lemma 7, there is to < oo large such that 

E (Ct(u) 1{0 < ( t (u) < F}) < | for all t > t . (13) 

To prove extinction of the frozen particles, we apply successively Proposition 2 to obtain the 
existence of an increasing sequence of times {tj : j > 0} tending to infinity such that 

P (Cs{u) = F for all s E (tj,t j+1 ) \ Q t . (u) = F) < - for all j > 0. 

In particular, since each time an F-site becomes an (F — l)-site, there are two particles that 
annihilate each other, the previous inequality implies that 

E(( tj (u)) - E(( tj+1 (u)) > 2x (l-^\ P(( tj (u) = F) = P(( tj (u) = F) (14) 

for all integers j > 0. To see this, note that if a site is frozen at time tj, then with probability at 
least 1/2 this site will be visited by an active particle before time tj+i, resulting in the annihilation 
of two particles. Therefore, the expected number of particles killed per site between times tj and 
ij+i is larger than the probability of a site being an F-site at time tj, which gives (14). Now, using 
(14) and then (13), we deduce that 

FxP(Ct J+ M = F) = E(( tj+1 (u)l{( tj+1 (u)=F}) < E(( t]+1 (u)) 

< E(C tj (u)) - P{Q tj {u) = F) 

= E(( t .(u)l{C tj (u) = F}) + E(C tj (u)l{C tj (u)^F}) - P(( t .( u ) = F) 

< FP(( t .(u) = F) + e/2 - P(( t .(u)=F) = (F - 1) P (( tj (u) = F) + e/2. 

In particular, a simple induction gives 

P(Ct„M = F) < (l-i) P(Ct.. 1 C)=F) + ^ 



i \ n -i / i \ n 

1 \ e 1 /. 1 \ e 



^ {'-f) + 2F 1-(1-1/F) = [ 1 -f) + 2 ^ 
for sufficiently large n. In view of Lemma 6, the previous inequality implies that 
lim P(( t (u) = F) = lim P (( t (u) = F) < e for all u G D. 

t— >oo n— >oo n 

Since this holds for all e > arbitrarily small, the lemma follows. □ 
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